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ABSTRACT 

We present the AKARI near-infrared (NIR; 2.5-5 /im) spectroscopic study of 
36 (ultra) luminous infrared galaxies [(U)LIRGs] at z = 0.01 — 0.4. We measure 
the NIR spectral features including the strengths of 3.3 /mi polycyclic aromatic 
hydrocarbon (PAH) emission and hydrogen recombination lines (Bra and Br/3), 
optical depths at 3.1 and 3.4 /im, and NIR continuum slope. These spectral 
features are used to identify optically elusive, buried AGN. We find that half 
of the (U)LIRGs optically classified as non-Seyferts show AGN signatures in 
their NIR spectra. Using a combined sample of (U)LIRGs with NIR spectra 
in the literature, we measure the contribution of buried AGN to the infrared 
luminosity from the SED- fitting to the IRAS photometry. The contribution of 
these buried AGN to the infrared luminosity is 5-10%, smaller than the typical 
AGN contribution of (U)LIRGs including Seyfert galaxies (10-40%). We show 
that NIR continuum slopes correlate well with WISE [3.4]— [4.6] colors, which 
would be useful for identifying a large number of buried AGN using the WISE 
data. 



Subject headings: galaxies: active - galaxies: ISM - galaxies: starburst - infrared: 
galaxies 
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INTRODUCTION 



Since the Infrared Astronomical Satellite (IRAS; iNeugebauer et al.l Il984j ) first opened 
the all-sky view of far-infrared universe, a large number of luminous infrared galaxies (LIRGs; 
10 11 ^ L IR (8-1000 /mi) < 10 12 L Q ) and ultraluminous infrared galaxies (ULIRGs; L m > 



10 12 LfT)) have been iden tified and studied extensively (see lSanders fc Mirabellll996l ; lLonsdale et al. 



2006; Soifer et al. 2008 for review 



In t he local Universe , many of them are interacting systems between gas- rich disk galax- 



ies | e.g 



2011 



Kim et al. 



2010a; but see 



2002 



Elbaz et al. 



Veilleux et al. 
20071 : 



2002 



Lotz et al 



Wang et al.lbood: iKavirailbood: iHwang et al. 



2008 



Ideue et al 



2009 



Kartaltepe et al 



2010 



for high-2 (U) LIRGs]. 
mass elliptical galaxies (e.g., 



2001 



Haan et al 



Tacconi et al.ll2002l ; iDasvra et al.ll2006 



They may evolve into quasars and then into intermediate- 



Sanders et al 



1988a 



Kormendv &; Sanders 



Veilleux et al.ll2009bl ; iRothberg fc Fischer 



1992 



Genzel et al. 



2010 



20111 ). Their enormous infrared luminosity comes from cool dust primarily 
heated by young stars [i.e., star formation (SF)] and hot dust heated by supermassive black 
holes (SMBHs) rapidly accreting matter [i.e., active galactic nuclei (AGN)]. Their contri 



bution to the infrared luminosity den sity increases with redshift (e.g.. lLe Floc'h et al.l 12005 



Magnelli et al. 


2009; 


Goto et al. 


2011) 



understand galaxy-galaxy interactions, starburst-AGN connection, and cosmic star forma- 
tion history. 

To identify the energy sources of galaxies (i.e., SF vs. AGN), the optical line ratios sensi- 
tive to the photoionizatiq n source have been used (e.g., 



Baldwin et al. 1981; Veilleux &; Qsterbrock 



19871 ; iKewley et al.ll2006[ ). Based on this me thod, the optical spectral types for a large num- 



ber of (U)LIRGs are also determined 



Goto 2005; Cao et al. 2006; Hou et al. 



e^. 



Veilleux et al 



2009 



Lee et al.l 



1995, 1999a 



Kewlev et all 12001 



201lh . However, the optical spectral 



classification can be uncertain, in particular for (U) LIRGs, due to the difficulty of detecting 
dust-enshrouded AGN. In this case, near-infrared (NIR) spectroscopy, less affected by dust 
extinction, is a more efficient tool. 

There are several ground-based f^-band (1.9-2.4 /mi) spectroscopic surveys searching for 
obscured AGN signatures in (U)LIRGs: the presence of broad Paa emission ce ntered at (rest- 



frame) 1.875 /mi or of high-excita t ion co ronal line [Si VI] at 1.963 /mi (e.g., I Veilleux et al. 



19971 . Il999ri Murphy et al.l Il999l |200l|). The ground-based L- band (2.8-4.2 urn) 



spec- 



troscopy was also used t o constrain their energy source (e.g., Ilmanishi fc Dudlevl Exx! 



Imanishi k Malonevl I2OO3I : Ilmanishi et~aD bood . boill ; iRisaliti et all bood . bold : ISani et al. 



20081 ). The SF-dominated (U)LIRGs show a strong emission line at 3.29 /im attributed 
to polycyclic aromatic hydrocarbons (PAHs), whereas AGN-dominated (U) LI RGs show a 
relatively PAH-free continuum attributed to larger-sized hot dust grains (e.g., iMoorwood 



-3 - 



19861 ; llmanishi fc Dudleyll2000l ). The strong absorption features at 3.05 and 3.4 /zm by H2O 
ice-covered dust grains and bare carbonaceous dust, respectively, are found in (U)LIRGs 
with buried AGN. However, these absorptions are weak or absent in normal SF (U)LIRGs 



where energy sources and dust are often spatially well mixed (e.g.. llmanishi fc Dudleyll2000 
Imanishi fc Maloneyil2003l ). For a similar reason, NIR continua of AGN-dominated (U) LIRGs 
can be much redder than those of SF-dominated (U)LIRGs (e.g., iRisaliti et al.ll2006l ). 



Thanks to the wide wavelength coverage (2.5-5 p,m ) and high sensitivity of the Infrared 



Cam era (IRC) on-board the AKARI space telescope ([Murakami et al.l l2007t lOnaka et al. 
20071). the L-band diagnos tic could be applied to more distant and faint galaxies (e.g. 



Imanishi et al.ll2008l . l2010af ). The AGN detection rate from L-band spectroscopy appears to 
be larger than that from Lf-band spectroscopy (roughly 70% vs. 20% in ULIRGs) because the 
iT-band diagnostic detects only obvious AGN. These NIR spectroscopic studies of (U)LIRGs 
suggest that there are many optically elusive buried AGN and that the buried AGN fraction 
increases with increasing infrared luminosity. 

In this study, we analyze the AKARI NIR spectra of 36 (U)LIRGd3 mainly from th e 
cross-correlation between the IRAS and Sloan Digital Sky Survey (SPSS; lYork et al.l l2000r) . 
Combining our new NIR spectroscopic data with those in llmanishi et al. r|2008l . boioah . we 
investigate the NIR properties of a large sample of (U)LIRGs. The structure of this paper is 
as follows. The target selection is explained in Section 2. Observations and data reduction 
are described in Section 3, and the method of NIR spectral analysis is given in Section 4. 
Our findings are discussed and summarized in Sections 5 and 6, respectively. Throughout 
this paper, we adopt flat ACDM cosmological parameters: H = 75 km s~ x Mpc -1 , Qm = 
0.3, and Q A = 0.7. 



TARGETS 



We obtained the AKARI NIR spectra of (U)LIRGs in three AKARI open-time pro- 
grams: the nature of new ULIRGs at intermediate redshift (NULIZ), NIR spectroscopy 
of composite and LINER LIRGs (CLNSL), and NIR spectroscopy of star- forming infrared 
galaxies (NISIG). The main ta r gets f or the NULIZ program were selected from a catalog of 
~320 ULIRGs in lHwang et al.l (120071). w hich was constructed by cross-correlating the IRAS 



faint source catalog (IMoshir et al. 



19921 ) with the galaxy redshift survey catalogs, including 



1 Although 5 out of 36 galaxies do not satisfy the definition of (U)LIRGs (i.e., Lir < 10 11 L & ), they are 
simply referred as (U)LIRGs in this study. We do not include them when we compare our results with those 
in previous studies. 
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the SDSS Data Re lease 4 (SPSS DR4; lAdelman-McCarthv et alJl2006f ). 2d F Galaxy Redshift 



Survey (2dFGRS; IColless et al.l 120011 ) . and 6dF Galaxy Survey (6dFGS; IJones et all 12004! . 
20051 ). These ULIRGs are detected up to z ~0.4. Because most ULIRGs observed with L- 



band spectrograph are at z < 0.2, we focus on 20 newly discovered ULIRGs at 0.2 < z < 0.4 
to extend the NIR spectroscopy to intermediate- z ULIRGs. The upper redshift limit ensures 
that the 3.3 /im PAH emission falls within the AKARI spe ctral coverag e . Add it ionally, four 



nearby ULIRGs were s electe d, as a control sample, fr om ISoifer et al.l (Il989l ). iLeech et al. 
Jl994h . lstanford et al.l J200oh . and lSanders et all J200J. 



For the CLNSL and NISIG programs, we first identified ~14,000 infrared galaxies by 
the cross-correlation of infrared sources in t he IRAS faint source c atal og with the spectro- 



scopic sample of galaxies in the SDSS DR7 (lAbazajian et al.l 120091 ; see iHwang et al.l 12010a 
for more details). We selected ~13,000 local infrared galaxies at 0.01 < z < 0.2. We 
then determined their optical spectral types [star-forming, (SF-AGN) composite, low ion- 
ization nuclear emissio n- line region (LINER ), an d Seyfert 2] usin g their emission line fluxes 
based on the criteria of Kewley et al. ( 20061 ) (see Lee et al. 2011 for more details). The line 
flux measurements were drawn from the Max-Planck-Institute for Astrophysics /Johns Hop- 
kins University value-added galaxy catalog (MPA/JHU VAGCs; Kauffmann et al. 2003a ; 
Tremonti et al. 2004 ; Brinchmann et al. 20041 ). We finally selected 60 i\~ s -bright§| non- Seyfert 
(i.e., SF, composite, and LINER) galaxies, covering a wide range of infrared luminosities 
(L IR = 10 10 -10 13 Lq). Higher priority was assigned to ^-bright galaxies in order to obtain 
high signal-to-noise ratio (S/N) AKARI spectra and to non-Seyfert galaxies in order to find 
"buried" AGN. Most targets are at high ecliptic latitudes due t o the limited sky visibility of 
AKARI satellite (Sun-synchronous orbit; iMurakami et al.ll2007l ). but this does not introduce 
any bias to the results. 

Among 48 observed targets, we use only 36 galaxies with median S/N per pixel for the 
continuum greater than three for the following analysis. Their basic information includ- 
ing the IRAS flux densities and optical spectral properties is summarized in Table [TJ We 
plot the infrar e d luminosity versus redshift of our sample in Figure [TJ The (U)LIRGs in 
Imanishi et al. ( 2008 , 2010a ). hereafter the Imanishi sampl^l, are also plotted for compari- 
son. The two samples are overall in a similar range of infrared luminosities, while our sample 



2 http://www. mpa-garching.mpg.de/SDSS/DR7/ 

3 -ftf s -band (2.17 fim) 20 mag arcsec -2 isophotal ellipse aperture magnitude drawn from the 2 micron 
all-sky survey (2MASS) < 14.5 mag 

4 This sample contains ULIRGs from the IRAS 1 Jy sample ( Kim fc Sanders! 1998 ) and L IRGs from the 



brigh t galaxy sample fBGS; ISoifer et al.lll987tlSanders et al.lll995l) and from the revised BGS (|Sanders et al 
20031) . 
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is more numerous than the Imanishi sample at 0.2 < z < 0.4. Imanishi et al. also preferen- 
tially selected non-Seyfert (U)LIRGs. In the result, non-Seyfert fractions in these samples 
are large compared to the fractions among their parent samples (from —55% to —70% in 
ULIRGs and from -65% to -80% in LIRGs). 



OBSERVATIONS AND DATA REDUCTION 



We observed the target galaxies with the IRC spectrograph (lOnaka et al.l 120071 ) on- 
board the AKARI satellite. It is a slitless spectroscopy, but the targets are located inside a 
lxl arcmin 2 window in order to avoid the spectral overlap with nearby sources. The 2.5-5 
fim wavelength range was covered at a spectral resolution of R —120 (—2500 km s _1 ) with 
NG grism. In the observing modes, the total on-source exposure time per pointing is about 
six minutes, and the 3a detection limit is roughly 0.6 mJy Each target was observed one 
to eight times (mostly twice), depending on its visibility and on i^-band magnitude. The 
observation log is given in Table HJ Note that the NULIZ data were obtained in Phase 2 
cooled by liquid helium, while the CLNSL and NISIG data were obtained in Phase 3 cooled 
by a cryocooler. 

Data reduction was carried out using the IDL packages prepared by the AKARI team: 
"IRC Spectroscopy Toolkit Version 20101025" for the NULIZ data and "IRC Spectroscopy 
Toolkit for Phase 3 data Version 20101025" for the CLNSL and NISIG data. Both pack- 
ages are available onlin^l This involves dark subtraction, linearity correction, flat-fielding, 
sub-frame combining, source detection, background subtraction and source extraction, wave- 
length calibration, and flux calibration. Because one pointing consists of eight (or nine) sub- 
frames, cosmic rays are efficiently removed. The wavelength calibration accuracy is —0.01 
/im (—1 pixel). The absolute flux calibr ation accuracy is —1 0% at the central wavelength 



and can be larger than 20% at the edge ( jOhyama et al.l 120071 ). 



We derived the spatial profile of each source in the two-dimensional spectra obtained 
with the slitless spectrograph. To extract one-dimensional spectra, we used 3-9 pixels (3-9 
x 1".46) along the spatial direction. The aperture width (see column 2 in Table 3) was 
determined to contain signals larger than the background noise level, which is the root- 
mean-square (rms) of fluxes at each side 11-15 pixel away from the center of the source. 
Since the aperture is large enough to cover the NIR emission even for extended sources, 
no aperture correction is applied. For the sources observed several times, the spectra were 
stacked by taking the median. The data collected during the late part of Phase 3 of AKARI 



5 |http://www.ir. isas.jaxa.jp/ASTRO-F/Qbservation/DataReduction/IRC7 
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are sometimes seriously affected by bad pixels and background noise because of the increase 
of the temperature in IRC. If there were some data sets substantially inferior to others, they 
were not used to obtain the final spectra. 

We display the final spectra of 36 (U)LIRGs with the continuum S/N > 3 in Figure 
|2j These spectra show various continuum shapes, ranging from steeply decreasing ones to 
nearly flat ones with increasing wavelength. In most spectra, the PAH emission feature at 
3.29 /im is clearly seen. Bra (2.63 /jm) and Bra (4.05 /im) hydrogen recombination lines 
and broad absorption features from H 2 ice (3.05 /xm) and bare carbonaceous dust (3.4 /im) 
are detected in some cases. 



4. MEASUREMENTS 



Figure [3] shows an example of spectral fitting to the NIR spectrum. We fit each emission 
line of the Bra, Br/3, and 3.3 /im PAH with a single Gaussian function after subtracting the 
local continuum defined by a linear fit around the line. Especially for the 3.3 /im PAH profile, 
we choose an asymmetric fitting range (3.15-3. 35 /im) in order to avoid contamination of 
the 3.4 /im sub-peak (see llmanishi et al.ll2010a| ). We compute the equivalent width of the 
PAH emission using the fitted Gaussian profile on top of the global continuum (see next 
paragraph). By using the global continuum instead of the local continuum, we can reduce 
the effect of adjacent absorption feature on the equivalent width measurement. We set an 
upper limit for the non-detection by assuming that the line has a Gaussian profile with 
a full width at half maximum (FWHM) and with an amplitude twice the local rms of 
the continuum-subtracted spectrum. As for the fixed FWHMs, we use the typical values 
measured from the AKARI spectra in this study: 2500 km s _1 for the Bracket lines and 
4000 km s -1 for the PAH feature. The linewidth of the Bracket lines is well matched to 
the instrumental resolution. The linewidth of the PAH feature is consistent with the PAH 



profile in other studies (e.g., type-1 sources in lTokunaga et al.lll99ll ). 



The global continuum slope V is determined using a single power-law model (F\ oc 
A r ) with data points less-affected by line features in the observed wavelengths at 2.7-4.8 
/im. The da ta points near the e dge of the spectra are excluded because of their large flux 
uncertainty (jOhyama et al.ll2007l ). The optical depths of the 3.1 and 3.4 /im dust absorption 
features (r 31 and r 3 4 , respectively) are calculated as the natural logarithmic ratios between 
the adopted continuum levels and the observed absorption profiles. The continuum levels are 
adopted from the global continuum, and the absorption profiles are obtained after smoothing 
the spectrum with a 10 pixel box (~0.1 /im) to minimize noise effect. We consider that the 
absorption feature is detected when the maximum difference between the continuum level 
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and the absorption profile is twice larger than the local rms. 

Measurements are listed in Table [31 The uncertainties estimated by considering the local 
rms values are also listed. The median S/Ns of detected features are 6.4, 5.2, 3.6, 3.3, and 
3.1 for 3.3 jttm PAH, Bra, Br/3 fluxes, T3.1, and r 3 . 4 , respectively. We present the equivalent 
width of 3.3 fj,m PAH emission line (EW3.3PAH) versus continuum slope (r) diagram in Figure 
HI It shows that most (U)LIRGs have EW 33 pah < 150 nm and — 3 < T < 0; the distribution 
of our sample is not different from that of the Imanishi sample. 



5. RESULTS AND DISCUSSION 



5.1. AGN Diagnostics 

5.1.1. AGN signature in the NIR spectra 



Emission at 3.3 /im is a prominent feature in (U)LIRG NIR spectra. It probably origi- 
nates from the reprocessing of ultraviolet (UV) radiation by PAH molecules. The contribu- 
tion of PfJ emission line at a similar wavelength is usually negligible. If (U)LIRGs host AGN, 
the EW3.3PAH is suppressed because hot dust emission from the AGN inc reases the NIR con - 
tinuum flux level and X-ray photons may destroy PAH molecules (e.g., ISmith et al.l 120071 ). 
A red NIR continuum (i.e., a high value of continuum slope) as well as strong absorption 
features at 3.1 and 3.4 /mi from H2O ice-covered dust grains inside molec ular clouds an d bare 
carbonaceous dust in the diffuse interstellar medium, respectively (see lDraindl2003l ). indi- 
cate the presence of highly obscured compact sources. While many of the se absorbed sources 



with red continuum slopes have been s hown to harbor buried AGN (e.g.. iRisaliti et al.ll2006 



spectral signatures (e.g.. iDesai et al 



Sani et al.ll2008l ; llmanishi et al.ll2010ah . deeply buried starbur s ts can also produce t he same 



2007 



Spoon et al.l 120071 ; IVeilleux et al.l l2009al ) and in 



general all that is required is a warm, highly obscured heating source. The small PAH 
equivalent width condition is useful for identifying weakly obscured AGN, while the large 
continuum slope and optical depth conditions are efficient in detecting highly obscured AGN. 
Therefore, all of these diagnostics are necessary to detect as many obscured AGN as possible. 



Following the criteria in llmanishi et al.l (j2010al ). we regard EW3.3PAH < 40 nm, V (F\ oc 



A r ) > —1, r 3 i > 0.3, and r 3A > 0.2 as AGN signatures, and classify sources satisfying at 
least one of these conditions as NIR AGN-detected galaxies (see column 10 in Table [3]). 
In the results, we find 19 AGN out of 36 (U)LIRGs based on these criteria. Among these 
AGN, there are 13 sources with EW3.3PAH < 40 nm (68%), five sources with T > — 1 (26%), 
five sources with r 3 .i > 0.3 (26%), and three sources with r 3 . 4 > 0.2 (16%). Note that 
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some sources satisfy more than one criterion. Therefore, EW3.3PAH < 40 nm is the primary 
criterion to select AGN. A similar trend is also seen in the Imanishi sample (56, 43, 45, and 
10% for EW3.3PAH, r, Tz.ii an d T3.4, respectively). 

We count the number of sources with AGN signatures among our sample (U)LIRGs in 
bins of optical spectral type and of infrared luminosity. These results are summarized in 
Table H] together with those from the Imanishi sample. The NIR AGN detection rates for 
our sample and the Imanishi sample are on average 53% and 51%, respectively, and agree 
well in each bin. Ilmanishi et al.l ( 120081 . |2010aj) found that the AGN signature in NIR spectra 
is more often detected in optical AGN-like and in more luminous infrared galaxies. Our 
sample appears to follow these trends, as shown in Figure [5j However, it is not conclusive 
with our data alone because of large uncertainties, in particular, for the NIR AGN detection 
rate depending on optical spectral type. In the combined sample of 180 (U)LIRGs with 
Lir, ^ 1O 11 L , the NIR AGN detection rate depends on optical spectral type as follows: 36% 
for SF (U)LIRGs, 55% for composite (U)LIRGs, and 66% for Seyfert 2 (U)LIRGs. Note that 
most of previously classif ied LINERs ar e calle d composites in this stu dy because we adopted 
the selection criteria of iKewlev et al.l (2006T) rather than those of IVeilleux fc Osterbrock 



(119871 ) that were used in Ilmanishi et al.l (120081 . l2010al ). There are two Seyfert 1 (U)LIRGs 
in the combined sample, and both of them show AGN signatures in their NIR spectra. The 
total NIR AGN detection rates for LIRGs and ULIRGs are 29% and 65%, respectively. If 
we select (U)LIRGs without any priority to non-Seyfert galaxies, the NIR AGN detection 
rates for (U)LIRGs increase slightly [e .g., 30% for LIRGs and 70% for ULIRGs based on the 



(U)LIRG sample in lYuan et al.l (I2010h ]. 



When we consider non-Seyfert galaxies with AGN signature in the NIR to be optically 
elusive buried AGN, the buried AGN fraction (i.e., the number ratio of buried AGN to non- 
Seyferts) increases with infrared luminosity: 24% for (U)LIRGs with L IR = 10 11 -10 12 L Q , 60% 
for (U)LIRGs with L m = 10 12 -10 12 - 3 L Q , and 84% for (U)LIRGs with L m = 10 12 - 3 -10 13 L Q . 
The higher AGN fraction in more luminous infrared galaxies has been reported in previous 
stud ies based on the data from not only AKARI but also the Infr ared Space Observatory (ISO; 



e.g 



mtz et al. 



1998 



Desai et al.l 



2007 



Rigopoulqu et al. 



Imanishi 



2009 



1999 



Tran et al. 



Valiante et al. 



2009; 



20011) and Svitzer Svace Telescope 



Veilleux et al. 2009a; Petric et al. 



201ll ). suggesting an important role of AGN in increasing infrared luminosity. 



As expected, we find many non-Seyfert (U) LIRGs with AGN signatures in their NIR 
spectra [buried AGN; 49% (55/113)]. On the other hand, a substantial number of Seyfert 
(U)LIRGs do not show the AGN signatures in the NIR [33% (14/43)]. The apertures used 
in the optical spectroscopy (our sample: 3" diameter fiber; Imanishi sample: 2" wide slit 
and 2 kpc extraction) are smaller than those of NIR spectroscopy (covers the entire galaxy 
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size). The line measur ements in the optica l spectra for our sample are aperture-corrected 
following the method in lHopkins et al.l ( 120031 ). but those for the Imanishi sample are not. The 
small aperture used for the optical spectral classification may miss the extended emission 
associated with star formation, and hence be more sensitive to weak, central AGN. This 
helps to understand the existence of Seyferts without AGN signatures in the NIR. We find 
no difference in the redshift distribution between Seyferts with and without AGN signatures 
in the NIR spectra. The different aperture size between the optical and NIR spectroscopy 
seems partially responsible for the disagreement of spectral types. 

Regardless of the aperture effect, if the hot dust emission from AGN is very weak 
because of a tiny covering of dust around the AGN, such AGN are clearly visible in the 
optical spectrum. However, they woul d be not distinguish able from SF-dominated galax- 
ies based on the NIR diagnostics (see iNardini et al.l I201CH ) . In contrast, when AGN are 
really heavily obscured, both optical and NIR diagnostics are less power ful. Then the ob 



servations at other wa velengths are nece ssary to detect them (e.g., X-ray: iBauer et al. 



2009a 



Teng fc Veilleuxll2010t mid /far- infrared: iFarrah et al.ll2007t ISpoon et al. 



Hatziminaoglou et al.ll2 010: El baz et al. 



2010b : radio: ISaiina et al.l 120081 : llmanishi et al.ll2010bj ). 



2011 



but see also Elbaz et al. 



2007 



2010: 



Veil 



2010 



eux et al. 



Hwang et al. 



In Figure O we present the optical AGN diagnostic diagram based on [Olll]A5007/H/3 
and [Nll]A6584/Ha line ratios, and NIR AGN diagnostic diagram based on EW3.3PAH and 
continuum slope. In panel (a), the optical AGN (Seyfert+LINER) without the NIR AGN 
signature have low [Olll]/H/3 ratios compared to those with AGN signatures both in optical 
and in NIR spectra. On the other hand, in panel (b), the NIR properties of optical AGN 
without the NIR AGN signature are not significantly different from those of non-AGN both 
in optical and in NIR spectra. 



5.1.2. AGN contribution to the infrared luminosity 

To measure the contribution of buried AGN to the infrared energy budget of (U)LIRGs, 
we use the infrared s pectral energy distribution (SED) templates and fitting routine of 



Mullanev et al.l ff201lh . DECOMPIEB. These templates consist of one AGN and five host- 
galaxy SEDs. For the host-galaxy SEDs, Spitzer mid-infrared spectra of starburst galaxies 
are extrapolated to the far-infrared using IRAS photometry. These host-galaxy SEDs are 
grouped into five categories, referred to as 'SB1' through 'SB5', in terms of their overall 
shape and relative strength of PAH features. For the AGN SED, the intrinsic SEDs of 



http: / /sites. google. com/site / decompir 
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AGN-dominated sources are derived after subtracting suitable host-galaxy components from 
the observed SEDs, and these SEDs are averaged. The infrared SEDs of AGN show a 
large spread, mainly dependent on dust distribution around AGN (i.e., smooth vs. clumpy 
torus structures). However, the different AGN SEDs do not significa ntly change the result- 



ing AGN contrib ution to the infrared luminosity in (U)LIRGs (see iMullaney et al.l 12011 



Pozzi et a 



20121). Based on th is SED fitting with sparse photometric data points such 



as IRAS, IMullaney et al.l (]201lh found that the intrinsic AGN luminosities measured are 
actually correlated with those from high-resolution mid-infrared observations of the AGN 
cores. 

We apply this routine to 69 (U)LIRGs with S/Ns > 3 at all four IRAS bands (12, 25, 
60, and 100 /im) in the combined sample, and fit their SEDs five times with AGN and one 
of host-galaxies by allowing renormalization of these two templates. We choose the best-fit 
solution with the lowest \ 2 value, computing the AGN contribution to the infrared (8-1000 
jum) luminosity from this template set. Figure [7] represents example SEDs with the best-fit 
AGN and host-galaxy templates. The AGN contribution in (U)LIRGs ranges from 0% to 
69%, and is on average 6-8% in LIRGs and 11-19% in ULIRGs. Because the combined 
sample preferentially includes non-Seyferts, the AGN co ntribution in this s tudy seems to be 

ULIRGs: -20% 



in 


Farrah et al. 


2003; 


in 


Risaliti et al. 


12010 



10% in Petric et al 



35-40% in lVeilleux et al.ll2009at -25% in lNardini et al.ll2010l : 15-20% 



2011 



Figure [8] shows the correlation of the AGN contribution with the presence of AGN 
signature in the NIR, optical spectral type, and infrared luminosity. Not surprisingly, the 
AGN contribution is higher in (U) LIRGs with AGN signature in the NIR than in those 
without AGN signat ure. The AGN con tribution clearly increases with increasing infrared 



luminosity (see also iDonoso et al.ll2012l ). similar to the trend of buried AGN fraction in 
Figure [5] (b). The AGN contribution for Seyfert (U)LIRGs is slightly larger than those for 
SF and composite (U) LIRGs. 



5.1.3. AGN diagnostics with WISE data 



Recently, the Wide- field Infrared Survey Explorer ( WISE: IWright et al.ll2010l ) opens up 
the opportunity to probe mid-infrared properties (3.4, 4.6, 12, and 22 fim) for a large sample 
of galaxies with excellent sensitivity. We use the WISE all-sky survey source catalogj to 
identify the WISE counterparts of the (U)LIRGs observed with AKARIIRC and (U)LIRGs 



7 |http : / / wise2 . ipac.caltech.edu /docs/ release / allsky / 
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nt 0.01 < z < 0.4 in the SDSS DR7 flHwang et alJl2010ah within 3" 



In Figure M (a ), we compare the continuum slope T with WISE [3.4]-[4.6] color (Vega 
magnitude system; iJarrett et al.ll201ll ) for the AKARI (U)LIRGs. We overplot the expected 
WISE [3.4]-[4.6] colors from simple power-law continuum models as a function of continuum 
slope (solid line). The continuum slope Y and WISE [3.4]-[4.6] color show a tight correlation 
(Spearman rank correlation coefficient = 0.82; the probability of obtaining the correlation 
by chance = 1.03 xl0~ 30 ) with some offset and scatter around the expected relation. Most 
sources have small values of continuum slope compared to the expectation because the pre- 
sented slopes (fitting range: 2.7-4.8 /xm) are affected by the spectrum at <3.4 /xm where 
the stellar population contribution is large (see lLee et al.ll2010f ). The scatter may come from 
the contamination by the 3.3 /iin PAH feature. For the galaxies at z < 0.13, the presence 
of 3.3 /iin PAH emission makes [3.4]-[4.6] color bluer than the color without the PAH emis- 
sion. On the other hand, if the galaxies are at 0.26 < z < 0.55, the presence of 3.3 jttm 
PAH emission makes [3.4]-[4.6] color redder than the color without the PAH emission. The 
amount of change in colors depends on redshift and on the strength of PAH emission. From 
the experiment with the spectra of our sample, we find that the presence of PAH emission 
can change the [3.4]-[4.6] color by ±0.2 mag, consistent with the scatter in Figure 9 (a). 

We find that the AGN selection criterion based on the continuum slo pe (i.e., T > —1) in 
this stud y is roughly equivalent to that of [3.4 ]— [4.6] > 0.8 suggested by IStern et al.l (120121 ) 



see also lAssef et al.l 12010 : IJarrett et al.l 1201 If ). If (U)LIRGs with AGN signatures in the 



NIR spectra are regarded as genuine AGN (filled symbols), the WISE color criterion selects 
AGN with 72% completeness (51 out of 71 genuine AGN satisfy the WISE color criterion) 
and 76% reliability (51 out of 67 objects which satisfy the WISE color criterion are genuine 
AGN). 

Figure IH] (b) shows the WISE [3.4]-[4.6] colors versus IRAS flux density ratios between 
25 and 60 fim (hereafter IRAS 25-60 /im colors) for the SDSS (U)LIRG sample. The IRAS 
25-6 Aim color is known to be associated with nuclei activity in infrared-luminous galaxies 



de Griip et allll985l : Zanders et aDll988bl ; iNeff fc Hutchingslll992l : IVeilleux et al 



2009a 



Lee et al.ll201ll ). AGN-dominated galaxies show warm IRAS 25-60 /im colors (/25//60 = 0.2) 
while SF-dominated galaxies show cool colors (/25//60 < 0.2). The composite and SF galaxies 
show similar distributions in this domain, but the AGN are significantly different. The 
median colors of AGN differ from those of non-AGN with significance levels of 6.8<r and 
7.4cx in the IRAS and WISE, respectively. Interestingly, there are a substantial number of 
SF (U)LIRGs in the lower-right corner. It seems real even if we consider large uncertainties 
associated with the IRAS colors. They have warm dust emission without hot dust emission, 
therefore appearing to be heavily obscured AGN. As a result, the WISE [3.4]-[4.6] color is 
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a good tracer of AGN-heated hot dust emission, but may not be sufficient to detect heavily 
obscured AGN. 



5.2. Comparison between Optical and Infrared Properties 

5.2.1. Star formation rate indicators 



The total infrared continuum, 3.3 fim PAH emission, and recombinati on lines includin g 
Ha and Bra of galaxies are useful indicators of star formation rate (SFR) ( IKennicuttlll998l ). 
In Figure [10] (a-b), we compare these SFR indicators: infrared luminosity versus (a) Ha 
and (b) 3.3 /im PAH l uminosities. The Ha luminosities are extinction-corrected using the 
Balmer decrement and ICalzetti et al.l (120001 ) extinction curve. The expected relationships 
are overplotted between these parameters (dott ed lines; hereafte r SF galaxy sequences): 
Lft a /Li R = 10~ 2 25 from the empirical relatio n in iKennicuttl (119981) and ^3.3pah/-^ir = 10~ 3 
from the observations of Mouri et al. ( 1990 ) and llmanishi ( 2002 ). These panels show that 
there are large offsets and scatters between the data and the expected relationships, even for 
SF (U)LIRGs without any AGN signature [pure SF (U)LIRGs; large filled symbols]. The 
offset from the SF galaxy sequence is larger in ULIRGs than in LIRGs, suggesting that Ha 
and 3.3 /im PAH emissions are more depressed in ULIRGs. 

There could be several reasons for the strong depression of Ha and 3.3 /im PAH emission 
in ULIRGs. (1) The amount of dust extinction could be systematically underestimated in 
ULIRGs. To check this effect, we plot the observed line ratios of Bra/Ha and Ha/H/3 in 
panel (c). The pure S F (U) LIRGs appear to have larger Bra/Ha ratios by considering 
the extinction curve of ICalzetti et al.l (120001 ) with a typical Ry value of 4.05. This can 
imply the need of a larger Ry value for the extinction correction in (U)LIRGs [e.g., Ry = 
5.01 denoted by the arrow]. Some studies also suggest the need of flatter /grayer extinction 
curves (i.e., large Ry valu es) in ULIRGs, which may be attribute d to sup ernovae- driven 



large -size dust grains (e.g.. iKawara et al.l 1201 ll : IShimizu et al.l 1201 ll : see also iBoquien et al. 



20121 ). However, because of large scatter in the data and different aperture size between 
the optical and NIR spectroscopy, it should be checked with a more extensive data set 
in future studies. (2) Even for ULIRGs without any AGN signature in optical and NIR 
observations, there could still be hidden AGN that play a role in the relative depression of line 
emission. To remove contamination of hidden AGN, it is necessary to search for AGN at other 
wavelengths. However, it is expected that the contribution of buried AGN is not significant 
as discussed in Section l5.1.21 (3) The line emission from ULIRGs may be intrinsically weak. 
Regardless of the presence of AGN, the star formation in ULIRGs produces large infrared 
continuum emission because a larger fraction of stellar UV photons is absorbed by dust 
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inside star- forming regions under the strong radiation field in ULlRGs ( lAbel et al.l 120091 ). 
The PAH emission in ULIRGs may be depressed in the sense that intense radiation fields do 
not produce photo-dissociation regions that are necessary for PAH emission, or l ead to the 



destruction of PAH carriers (see IVoitlll992l ; ISmith et al.ll2007l : IVeilleux et al.ll2009af ) . It is still 



unclear whether t he ionization state of grains is actually related to infrared luminosity (e.g., 



Desai et al.ll2007t ISmith et al.ll2007t IVeilleux et al.ll2009at llmanishi et al.ll2010al : IPetric et al. 
201lh. 



5.2.2. Dust extinction as a function of infrared color 



In Figure [HI we present the Bra/H/3 line ratio of (U)LIRGs as a function of IRAS 25- 
60 fim color. There is an anti-correlation between the two parameters with Spearman rank 
correlation coefficient = —0.67 and the probability of obtaining the correlation by chance 
= 0.02. This suppor t s that warmer sources are less extinguished t han cooler sources (see 



Veilleux et al.l Il999al l2009al : Keel et al.l 120051 : llmanishi et al.l 120081 ). When using Balmer 



decrement values instead, such a correlation becomes weaker since the optical lines do not 
trace well the buried sources. 



6. SUMMARY 

We obtained AKARI 2.5-5 fim spectra of 36 (U)LIRGs, selected mainly from the IRAS- 
detected galaxies in the SDSS. We measured the NIR spectral features including continuum 
slope, 3.3 fim PAH strength, optical depths at 3.1 and 3.4 fim, and Bracket lines in order 
to find AGN signatures. Together with samples in the literature and ancillary data, we 
compared optical and infrared properties of (U)LIRGs. Our primary results are summarized 
below. 



1. We found that 52% (14/27) of optical non-Seyfert galaxies in our (U)LIRG sample 
show the AGN signature in their NIR spectra. The NIR AGN detection rate for the 
combined sample is higher in composite (U)LIRGs than in SF (U)LIRGs, and increases 
with infrared luminosity. 

2. We fitted the IRAS photometric data for 69 (U)LIRGs with the AGN/starburst SED 
templates to compute the AGN contribution to the infrared luminosity. We found that 
the contribution of buried AGN to the infrared luminosity is 5-10%, smaller than the 
typical AGN contribution in (U)LIRGs including Seyfert galaxies (10-40%). 
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3. The NIR continuum slopes correlate well with WISE [3.4]-[4.6] colors. Using the WISE 
[3.4]-[4.6] color vs. IRAS 25-60 /iin color domain, we found sources associated with 
warm dust emission, but without hot dust emission. The WISE color is useful for 
identifying a large number of AGN, while it can miss heavily obscured AGN. 
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Fig. 1. — Infrared luminosity vs. r edshift distribut i on of (U)LI RGs observed with AK ARI 
IRC. The (U)LIRGs in this study, llmanishi et al.l ( 120081 ) . and llmanishi et al.l ( ]2010al ) are 
represented by black filled circles, red open triangles, and blue open squares, respectively. 
In the upper and right panels, the redshift and luminosity distributions of (U)LIRGs are 
shown, respectively (our sample: solid line with filled circles; Imanishi sample: dotted line 
with open diamonds). 
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Fig. 2. — AKARI NIR spectra of 36 (U)LIRGs in our sample. For each object, their optical 
spectral type (SF=star-forming; Co=composite; LI=LINER; S2=Seyfert 2; Sl=Seyfert 1; 
Un=Unclassified) and redshift are shown. The spectra presented here are smoothed using a 
3 pixel box (~0.03 /im) and are normalized to peak values, denoted in units of 10~ 13 ergs s~ x 
cm -2 /im -1 . The dotted lines indicate the wavelengths of Br/3 (2.63 /im), H 2 ice-covered 
dust absorption (3.05 /im), PAH emission (3.29 /im), bare carbonaceous dust absorption (3.4 
/im), and Bra (4.05 /im) features. 
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Fig. [2j — Continued 




Fig. 3. — Example of spectral fitting for F10565+2448. The background solid line is the 
observed spectrum. The Gaussian fit lines to measure 2.63 fim Br/3, 3.3 /xm PAH, and 4.05 
fim Bra emissions are overplotted. The solid-dotted power-law line represents the global 
continuum. The wavelength intervals not used in the continuum fit are denoted by dotted 
lines. The thick line around 3.1 //m is the 10 pixel smoothed profile of 3.1 /im absorption. 
The arrow shows the wavelength where the difference between the continuum level and 
absorption profile is maximized. In this case, 3.4 /im absorption is not detected. 
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Fig. 4.— EW3.3PA H vs. continuum slope T d iagram. The (U)LIRGs in this study, 



Imanishi et al.l (120081 ) . and llmanishi et al.l ( j2010af ) are represented by black filled circles, 



red open triangles, and blue open squares, respectively. The sources with the upper limits of 
EW3.3PAH are represented by arrows. The continuum slope T was defined by F x oc A r in this 
study, while it was defined by F v oc A r in Imanishi et al. Thus, r t hi S study equals ri manis hi — 2, 
and the values of r tn is study and ri manisn i — 2 are presented for fair comparison. The EW3.3PAH 
and continuum slope distributions are shown in the upper and right panels, respectively (our 
sample: solid line with filled circles; Imanishi sample: dotted line with open diamonds). In 
the upper panel, the (U)LIRGs with EW3.3PAH upper limits are not included. 
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Fig. 5. — (a) AKARI-based AGN detection rate as a function of optical spectral type (See 
Section 15.1.11 for the definition of AGN signatures in the NIR spectra), (b) Buried AGN 
fraction as a function of infrared luminosity. The buried AGN fraction means the number 
fraction of NIR AGN among optical non-Seyfert galaxies. The blue circles, red diamonds, 
and black stars indicate the results from our sample, Imanishi sample, a nd the combined 



samp le, respectively. The error bars are based on Poisson statistics (see iDe Propris et al. 
2004h . 
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Fig. 6. — (a) Optical diagnostic dia gram based on [ Qui] A 5007/ H/? and [Nll|A6 5 84/Ha line 
ratios. The (U)LIRGs in this study, llmanishi et al.l (120081 ) . and llmanishi et al.l (j2010al ) are 
represented by circles, triangles , and squares, respectively. Th e optical line r atios of the 
Imanishi sample are taken from I Veilleux et al.l ( 119951 . Il999af ) and iKewley et al.l (j200lf ). The 
sources with (/without) AGN signatures in their NIR spectra are denoted by re d filled (/blue 



open ) symbols. The solid and dashed lines indicate the maximum starburst (IKewley et al. 



20011 ) and pure star formation (iKauffmann et al.ll2003bl ) lines, respectively, (b) NIR diag- 



nostic diagram based on EW 3 . 3 pah and continuum slope T. The optical SF, composite, and 
AGN (LINER + Seyfert) sources are denoted by small blue open, large green open, and large 
red filled symbols, respectively. The (U)LIRGs optically unclassified or with upper limits are 
not included. The dotted line shows the AGN selection criteria in this study: EW 3 . 3 pah < 
40 nm or r > — 1. 
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Fig. 7.- The best fit SEDs for (a) F1813 1+6820 and (b) F134 51+1232 derived using 
the infrared SED-fitting routine, DECOMPIR jMullanev et al.lboilh . The filled circles with 
errorbars indicate their four IRAS flux densities. The black solid, blue dotted, and red 
dashed lines represent the total, host-galaxy, and AGN SEDs, respectively. The labels in the 
upper left indicate which host-galaxy template is used and the derived AGN contribution to 
the total infrared luminosity. 
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Fig. 8. — The AGN contribution to the infrared luminosities of (U)LIRGs in bins of (a) 
optical spectral type, (b) presence of AGN signature in the NIR spectra, and (c) infrared 
luminosity. The mean values with sampling errors are shown. 
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Fig. 9.— (a) Comparison between the WISE [3.4]-[4.6] color (Vega) and NIR contin- 
uum slope T for (U )LIR GS observed wi t h AK ARI IRC. The (U)LIRGs in this study, 
Imanishi et al.l (120081 ) , and llmanishi et al. ( boiOah are represented by circles, triangles, and 



squares, respectively. The optical composite and AGN (/SF) sources are denoted by large 
(/small) symbols, and the sources with (/without) AGN signature in NIR spectra are shown 
by filled (/open) symbols. The symbols are also color-coded according to optical spectral 
types (blue: SF; green: composite; red: AGN). The solid line represents the WISE colors 
expected from simple power-law continuum models (F\ oc A r ). The vertical dotted and 
horizontal dash ed lines indi c ate th e criteria separating AGN from SF-dominated galaxies in 
this study and Istern etaJ (boij ). respectively, (b) The WISE [3.4]-[4.6] color vs. IRAS 
flux d ensity ratio of 25 to 60 /xm diagram. The small open symbols are (U)LIRGs in SDSS 
DR7 (IHwang et al.ll2010al ). Their median values are denoted by large filled symbols with 
sampling errors. The red circles, green squares, and blue triangles are for AGN, compos- 
ite and SF galaxies determined in the optical spectra. The (U)LIRGs with S/Ns > 3 at 
(WISE) 3.4, 4.6, (IRAS) 25, 60 /im are only presented. The error bars in the upper-left 
corner indicate the typical uncertainties associated with the colors. The vertical dotted and 
horizontal dashed line s ind icate the cri t eria s eparating AGN from SF-dominated galaxies in 



Sanders et al.l (jl988bl ) and lStern et al.l (120121 ). respectively. 



31 



13 



hJ 12 

(J 

? 11 



10 



1 1 1 1 1 


1 I 1 1 1 

-(b) ; 






■ v . ;r . 


■ v 






• 

o 

■ ■ ■ ■ ' ■ ' 


• 

■ ■ : . . . i i i 


7 8 9 10 


6 7 8 9 


!og (L H a/L©) 


lo g (L 3 .3pah/L©) 



8 

o 



1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 



1 1 1 1 1 1 

■ (c) 


■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

R„=4.05 ■ 
R v =5.01 1 

i\ A. 




.□ ■ . / .. c j 






□ 


: : ''o" 


; ♦ 

. i . . . i 





0.4 0.6 0.8 1.0 1.2 1.4 1.6 

lOg (W f H P ) 



Fig. 10. — Comparisons of infrared luminosity with (a) H a and (b) 3.3 <im PAH lumi nosi- 
ties. The (U)LIRGs in this study, llmanishi et all ( 120081 ) . and llmanishi et all ( ]2010ai ) are 
represented by circles, triangles, and squares, respectively. The optical SF (/composite and 
AGN) sources are denoted by large (/small) symbols, and the sources without (/with) AGN 
signature in NIR spectra are shown by filled (/open) symbols. The symbols are also color- 
coded according to optical spectral types (blue: SF; green: composite; red: AGN). In each 
panel, the dotted line indicates the SF galaxy sequence (see Section I5.2.ip . (c) Comparison 
between the observed Bra/Ha and Ha/H/3 line ratios. The symbols are the same as in 



panels (a-b), but the diamond in the 
galaxies in this plane (c ase B theory; 
of Ay = 8 mag from the ICalzetti et al. 



ower-left corner indicates th e intrinsic position of SF 



Osterbrock fc Ferlandll2006l ). The reddening vectors 
(boooh extinction curve with Ry = 4.05 and 5.01 are 
shown by dotted line arrows. Since the curve is given at shorter wavelengths than 2.2 /zm, 
the e xtinction at 4.05 jum (A-Rm) is extr apolated fro m A 2 .2/ ; ,m using the mod el of A\ cc A~ L75 
(e.g., ICardelli et al.lll989t lDraindll989l : but see also iNishiyama et al.ll2009l ). We fit a linear 



relation to the data of pure SF (U)LIRGs after removing three outliers with 2a clipping 
(open star symbols), and so obtain the curve with Ry = 5.01. 
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Fig. 11. — The Bra/Ha line ratio vs. IRAS ftux density ratio between 25 and 60 /zm for the 
AKARI (U)LIRG sample. The symbols are the same as in Figure 9 (a). 



Table 1. Basic information and optical spectral properties. 
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.510 


12.53 


Co 


5.12 


9. 


34 





.95 





.59 


1/2 


F16533+6216 


16 


53 


52. 


1 


+62 


11 


50 


0. 


106 


<0. 


.059 


0. 


169 


1. 


480 


2 


.503 


11.94 


Co 


5.75 


9. 


34 





47 





.52 


1/2 


F22509-0040 


22 


53 


33. 


.0 


-00 


24 


43 





058 


<0 


174 


0. 


716 


5. 


143 


5. 


028 


11.87 


Co 


6.34 


9. 


07 





.66 





.63 


1/2 



CO 
CO 



Table 1 — Continued 



Object 


RA 


Dec 


z 


/12 




/(in 


/ioo 


log Lir 


Type 


Ha/H/3 


log L Ho 


[Om]/H/3 


[Nil] /Ha 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


F23223+1459 


23 24 49.4 


+ 15 16 32 


0.014 


0.144 


<0.339 


1.310 


4.461 


10.29 


LI 


3.89 


7.13 


1.56 


1.46 


1/2 


09022-3615 


09 04 12.7 


-36 27 01 


0.060 


0.211 


1.154 


11.470 


11.430 


12.22 


SF 


5.67 


9.27 


1.16 


0.33 


12/4 



Note. — Column 1: Object name in the IRAS catalogs. Columns 2-3: Right ascension and declination in units of ™ m s and ° ' ", respectively (J2000). Column 4: 
Redshift. Columns 5-8: The IRAS flux densities at 12, 25, 60, and 100 (im (Jy). Column 9: Logar ithm of infrared l umino sity (Lq), derived from the IRAS fluxes using 
the formula, Ltr = 2.1 X 10 39 X D(Mpc) 2 X (13.4 8 X /i 2 + 5.16 X / 25 + 2.58 X / 60 + /ioo), in lSanders fc Mirabel jl996h . For sources with upper limits, we follow the 
method described in llmanishi et alj (I20081. l2010ah . The upper and lower limits on the infrared luminosity are obtained by assuming that the actual flux is equal to the 
IRAS upper limit and zero value, respectively, and the average of these values is adopted as the in frared luminosity. Co lumn 10: Optical spectral type, determined with 
spectroscopic data covering at least from H/3+[Om] A5007 to Ha+[Nn]A6584 using the criteria of iKewlev et al. | l l2006h (SF=star-forming; Co=composite; LI=LINER; 
S2=Seyfert 2; Sl=Seyfert 1; Un=Unclassi fied). Column 11: O bserved Ha/H/3 ratio. Column 12: Logarithm of Ha luminosity (Lq), corrected for dust extinction 
using the observed Ha/H/3 ratio with the [C^lzeJiti^eJ^^l. j|2000|) extinction law (Ry = 4. 05) on assum i ng the intrinsic Ha/H/3 ratio of 2.85 for SF galaxies and 3.1 
for composite plus A GN galaxies JOsterbrock fe Ferlandl l2006h . For SDSS (U)LIRGs in iHwang et alj feOlOal l, the aperture correction is also applied following the 
method suggested by [Hopkins et al .1 l|2003l l based on the difference between fiber and Petrosian magnitudes to minimize the small fixed-size (3") aperture effect, if 
available. Column 13: [Oni] A5007/H/3 r atio, corrected for du st extinction. Colu mn 14: [Nll1A 6584/Ha r atio, corrected for dust extinction. Col umn 15: References of ht^ 
the identification/optic al p roperties. (lllHwang et al.ll2010al; (21 T his st udy: (3llHwanget aLll2007l: (41 iLee et alj|201ll; (5 1 IStanford et al.ll2000l: (6) ISoifer et af]|l989l; 
(7l lVeilleux et al.lll995l: f8) lYuan et"al]|2010l: (91 iLeech et al.lll994l: (101 IVeron-Cettv fc Veron II2OO6I: (lll lHou et al.ll2009l : (12l lSanders et aljbooa . 
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Table 2. Observation log. 



Object 


Program 


Observation ID 


Observation Date 


F01159-0029 


CLNSL 


3350035-001- 


-002 a 


2009-07-09 


F01212-5025 


NULIZ 


3050012-001- 


-002 


2006-12-12 


F02595-4714 


NULIZ 


3050008-001- 


-002 


2007-01-07-08 


F03130-3119 


NULIZ 


3050009-001- 


-002 


2007-01-25 


F03483-4704 


NULIZ 


3050010-001- 


-002 


2007-01-22 


F08082+2521 


NISIG 


3610001-001 




2009-10-22 


F09019+5136 


NISIG 


3610011-001- 


-002 


2009-10-25-26 


F09174+0821 


NISIG 


3610035-001- 


-002 


2009-11-12-13 


F09395+3939 


CLNSL 


3350003-001- 


-003 


2009-05-05 


F09501+5535 


NULIZ 


3050001-001- 


-002 


2006-11-01, 2007-04-29 


F10035+4852 


NISIG 


3610008-001- 


-002 


2009-11-06-07 


F10300+3509 


NULIZ 


3051003-001- 


-002 


2007-05-17 


F10565+2448 


NULIZ 


3051019-001 




2007-05-28 


F11087+1036 


CLNSL 


3350028-001 




2009-06-05 


F11163+5707 


NISIG 


3610013-001- 


-004 


2009-11-14 


F11213+6556 


NULIZ 


OUUUUOU UU -L 




2007-05-04 


F11514+1212 


NULIZ 


3051008-001- 


-002 


2007-06-15 


F12207+6329 


CLNSL 


3350024-001- 


-002 


2009-05-14 


F12232+5532 


NULIZ 


3051006-001- 


-002 


2007-05-25 


F12469+4359 


NULIZ 


3051010-001- 


-002 


2007-06-11 


F12471+4759 


NULIZ 


3051009-001- 


-002 


2007-06-07-08 


F12479+3705 


NULIZ 


3051005-001- 


-002 


2007-06-16 


F12536+5113 


NISIG 


3610028-001- 


-004 


2009-12-06-07 


F13099+4627 


CLNSL 


3350015-001 




2009-06-12 


F13300+6652 


CLNSL 


3350025-001- 


-008 b 


2008-11-19-22, 2009-05-18-20 


F13380+3339 


NULIZ 


3051001-001- 


-002 


2007-06-30 


F14014+3718 


NULIZ 


3051007-001- 


-002 


2007-07-03 


F14166+6514 


NULIZ 


3051014-001- 


-002 


2007-05-28 


F14379+5420 


NULIZ 


3050004-001- 


-002 


2006-12-24 


F14413+3730 


NULIZ 


3050005-001- 


-002 


2007-01-11-12 


F14541+4906 


NULIZ 


3050006-001- 


-002 


2007-01-04 


F15002+4945 


NULIZ 


3051015-001- 


-002 


2007-07-06 


F16533+6216 


CLNSL 


3350023-001- 


-002 


2009-01-08 


F22509-0040 


CLNSL 


3350036-001- 


-002 


2009-06-04 


F23223+1459 


CLNSL 


3350037-001- 


-005 


2009-06-19-20 


09022-3615 


NULIZ 


3051018-001- 


-002 


2007-05-26 



a The 2nd pointing was not used to obtain the final spectra. 
b The 4th pointing was not used to obtain the final spectra. 
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Table 3. Measurements and AGN signature. 



Object 


Aper. 


h 


.3PAH 




EW3.3PAH 


/ 


Bra 


/Br/3 




T3.1 






T3.4 






r 






AGN 


(1) 


(2) 




(3) 






(4) 


(5) 










(7) 






(8) 






(9) 






(10) 


F01159-0029 


8.8 


6 


Q 

y 


_j_ 
33 


1 


8 


59 ± 16 




< 


2.0 


2.6 ± 


1 


1 




< 





10 






-2 


70 
1 Z 


_l_ 
± 





07 


X 


F01212-5025 


10.2 






< 





5 


< 13 








< 





5 









15 


< 


0.14 


-2 


uo 


_1_ 
AZ 





06 


O 


F02595-4714 


8.8 





o 
o 


33 





1 


18 ± 6 








0.8 ± 





1 









25 


< 


0.35 


-2 


1 ^ 


-L 
AZ 





07 


O 


F03130-3119 


8.8 


1 


9 
Z 


_|_ 
33 





1 


81 ± 9 








< 





6 









39 


< 


0.15 


-2 


00 


_|_ 

33 





10 


X 


F03483-4704 


8.8 






< 





5 


< 28 








< 





6 




< 





43 


< 


0.24 


-1 


00 


_l_ 
33 





10 





F08082+2521 


13.1 


6 


a 
D 


_|_ 

33 


1 


7 


40 ± 10 





















13 






-3 


"1 "1 


_1_ 
33 





08 


X 


F09019+5136 


10.2 


15 


4 


33 


2 


4 


104 ± 16 






1.8 


< 


1 


4 










12 






-2 


^7 


-1- 

33 





06 


X 


F09174+0821 


11.7 


10 


i 

i 


_|_ 

33 


1 


6 


48 ± 7 






1.0 


< 


2 


1 




<■ 





10 






-2 


00 


33 





05 


X 


F09395+3939 


7.3 


1 


9 


_|_ 

33 





2 


60 ± 8 






0.7 


< 





8 




<c 





13 






-2 


lo 


_|_ 

33 





14 


X 


F09501+5535 


7.3 





r: 
o 


33 





1 


58 ± 8 








< 





4 


u.yu 


-L- 





16 


0.59 ± 


0.30 


-1 


on 
uu 


3_ 
33 





09 





F10035+4852 


8.8 


14 


a 
D 


_|_ 
33 


1 


3 


83 ± 8 






1.3 


< 


2 












13 






-2 


■. ) 


_|_ 

33 





05 


X 


F10300+3509 


7.3 





c 

u 


33 





1 


43 ± 11 








< 





6 




< 





38 






-2 


01 


_L 
33 





15 


X 


F10565+2448 


10.2 


37 


n 
\) 


33 


5 


6 


111 ± 17 


A A 


33 


0.8 


4.9 ± 


1 


4 


U.OO 







06 






-2 


ou 


_1_ 
33 





05 





F11087+1036 


10.2 


5 


D 


_|_ 
33 


1 


1 


88 ± 18 






0.8 


< 


1 


6 









14 


< 


0.14 


-3 


U / 


_|_ 

33 





11 


X 


F11163+5707 


7.3 


1 


r; 
□ 


-1- 

31 





3 


66 ± 14 






0.9 


< 


1 













15 






-1 


Q7 


3_ 
33 





15 


X 


F11213+6556 


7.3 











5 


< 37 








< 





8 




<■ 





26 


< 


0.27 


-2 


zo 


_|_ 

33 





16 





F11514+1212 


5.8 











5 


< 61 








< 





5 









15 






-A 


OA 


_|_ 

33 





22 


X 


F12207+6329 


10.2 


15 


7 


± 


1 


8 


180 ± 20 


3.0 


± 


0.3 


2.3 ± 





4 




< 





09 






-2 


79 


± 





08 


X 


F12232+5532 


8.8 


1 





± 





1 


39 ± 3 








< 





7 




< 





16 


< 


0.21 


-1 


86 


± 





08 





F12469+4359 


7.3 






< 





6 


< 88 








< 





5 


1.21 


± 





47 






-3 


45 


± 





19 





PI 9471 -\-d7 r *Q 


7.3 





6 


± 





2 


22 i 7 








< 





5 




< 





20 




0.14 


_2 


01 


± 





07 




F12479+3705 


7.3 





6 


± 





2 


28 ± 8 








< 





6 


0.22 


± 





07 


< 


0.11 


-2 


53 


± 





10 





F12536+5113 


7.3 


4 





± 





6 


147 ± 23 




< 


1.6 


< 


1 


2 




< 





25 






-2 


62 


± 





28 


X 


F13099+4627 


11.7 






< 


5 


1 


< 15 




< 


1.2 


< 


3 


2 


0.12 


± 





04 






-3 


20 


± 





05 





F13300+6652 


7.3 






< 





3 


< 22 




< 


0.3 


< 





5 




< 





24 


< 


0.23 


-2 


93 


± 





11 





F13380+3339 


7.3 


1 


1 


± 





1 


44 ± 5 








< 





6 




< 





35 


< 


0.15 


-2 


16 


± 





08 


X 


F14014+3718 


7.3 


1 


4 


± 





1 


76 ± 5 








< 





7 




< 





14 






-1 


68 


± 





08 


X 


F14166+6514 


5.8 






< 





7 


< 31 








< 





5 




< 





18 






-2 


04 


± 





20 





F14379+5420 


4.4 





6 


± 





2 


71 ± 27 








< 





3 




< 





42 






-1 


86 


± 





21 


X 


F14413+3730 


7.3 





9 


± 





3 


55 ± 20 








< 





5 




< 





37 


< 


0.12 


-2 


54 


± 





Id 


X 


F14541+4906 


8.8 






< 





4 


< 3 








< 





7 




< 





08 


0.29 ± 


0.05 





62 


± 





07 





F15002+4945 


7.3 






< 





7 


< 3 








< 





7 




< 





11 






-0 


52 


± 





05 





F16533+6216 


10.2 


6 


7 


± 





8 


137 ± 16 


0.6 


± 


0.3 


0.4 ± 





2 


0.42 


± 





14 






-2 


33 


± 





08 





F22509-0040 


7.3 


17 


8 


± 


1 


9 


125 ± 15 


3.0 


± 


1.2 


< 


2 


1 


0.62 


± 





19 


0.38 ± 


0.14 


-0 


91 


± 





08 





F23223+1459 


13.1 






< 


1 


5 


< 3 




< 


0.9 








0.23 


± 





06 






-3 


38 


± 





05 





09022-3615 


13.1 


21 





± 


4 


4 


44 ± 10 


4.9 


± 


0.9 


4.3 ± 


1 


2 




< 





08 


< 


0.06 





55 


± 





12 






Note. — Column 1: Object name in the IRAS catalogs. Column 2: Aperture width used for the extraction ("). Column 3: Observed 
flux of the 3.3 /im PAH emission (10~ 14 ergs s _1 cm -2 ). Column 4: Rest-frame equivalent width of the 3.3 /im PAH emission (nm). 
Columns 5-6: Observed fluxes of Bra and Br/3 (10 -14 ergs s" 1 cm -2 ). The Bra and Br/3 fluxes of sources at z > 0.2 and z < 0.02 
are not presented, respectively, because the measurements near the edge of the spectra are uncertain. Columns 7-8: Optical depths of 
the 3.1 fim H2O ice and 3.4 (im bare carbonaceous dust absorption features. The 3.4 /im optical depths are not presented if they are 
significantly affected by the 3.4 (im PAH sub-peak. Column 9: Continuum slope T (F\ oc A r ). Column 10: AGN signature from the 
NIR features (EW 3 . 3 pah < 40 nm, T3.1 > 0.3, T3.4 > 0.2, or F > -1). 



Table 4. NIR AGN detection rate in bins of optical spectral type and of infrared 

luminosity. 



Sample 


log Lir 






Optical spectral type 








(Lq) 


Star-forming 


Composite 


LINER 


Scyfert 2 


Seyfert 1 


Unclassified 


This study 


12.3-13.0 
12.0-12.3 
11.0-12.0 
10.0-11.0 


100% (1/1) 
50% (2/4) 
0% (0/3) 
0% (0/3) 


100% (3/3) 
60% (3/5) 
50% (3/6) 
... (0/0) 


... (0/0) 
... (0/0) 
... (0/0) 
100% (2/2) 


50% (1/2) 
... (0/0) 
100% (1/1) 
... (0/0) 


100% (1/1) 
... (0/0) 
... (0/0) 
... (0/0) 


0% (0/1) 
50% (2/4) 
... (0/0) 
... (0/0) 


Imanishi+08,10 


12.3-13.0 
12.0-12.3 
11.0-12.0 


100% (2/2) 
67% (4/6) 
18% (3/17) 


77% (10/13) 
61% (19/31) 
25% (5/20) 


... (0/0) 
50% (1/2) 
... (0/0) 


75% (6/8) 
70% (14/20) 
50% (5/10) 


... (0/0) 
... (0/0) 
100% (1/1) 


0% (0/2) 
57% (4/7) 
20% (2/10) 



Note. — The values in parenthe ses mean the number of sour ces with AGN signature/total number of sources in each bin 
from our sample a nd th e sample in llmanishi et all {2008, 2010 a|), The opti cal spectral types of the Imanishi sample are taken 
from I Yuan et al.1 feoich so that the same criteria of iKewIe^etal] i2006l ) are used for the two samples. For sources with 
multiple nuclei, if AGN signatures are seen at least in one nucleus, then these sources are considered to be AGN. For an 
unbiased comparison, the additional interesting sources in the Imanishi sample are not included. 



